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ABSTRACT 

A special form of the isotropic metric in cylindrical coordinates is used to construct 
what may be interpreted as the General Relativistic versions of some well known 
potential-density pairs used in Newtonian gravity to model three-dimensional distri- 
butions of matter in galaxies. The components of the energy- momentum tensor are 
calculated for the first two Miyamoto-Nagai potentials and a particular potential due 
to Satoh. The three potentials yield distributions of matter in which all tensions are 
pressures and all energy conditions are satisfied for certain ranges of the free param- 
eters. A few non-planar geodesic orbits are computed for one of the potentials and 
compared with the Newtonian case. Rotation is also incorporated to the models and 
the effects of the source rotation on the rotation profile are calculated as first order 
corrections by using an approximate form of the Kerr metric in isotropic coordinates. 
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1 INTRODUCTION 

Analytical axially symmetric solutions of gravitating matter provide very accurate models for the mass distribution of galax- 
ies. In Newtonian theory, models for globular clusters and spherical galaxies were presented by iPlummer. QQllh an djKinel 
il966l) . On the other hand, to describe the mass distribution within highly flatt ened axisymme tric galaxies. IXoomr j Im63) 
foun d a family of den sity potentia l pairs, the first of th em was flrst derived bv iKuzminI il956l) . Later, Miyamoto and Na- 
gai ijMivamoto fc Nagai, .1975: Nagai fc Mivamotol llfl?^ ) "thickened- up" Toomre's series of disk models to yield pairs of 
three-dimensional potential and density functions. In a similar way, ISatohl (|l98(t obtained a family of three-dimensional 
axisymmetric mass distributions from the higher order Plumme r models. For a detailed account and the description of other 
potential-density pairs used in galactic models, see for example iBinnev fc Tromaina 1.198 7'). 

In General Relativity, several exact d isk-l ike solutions have been fou nd, among them the static disks without radial pres- 
sure studied bv lBonnor fc Sackfieldl lll96^ and lMorgan fc MorganI lll969l^ , and disks with radial pressure bv lMorgan fc MorganI 
Several classes of exact solutions of the Einstein fi eld equations corresponding to static thin disks with or with- 
out r adial pressure have been obtained by different autho rs llLvnden-BeU fc Pineaullll978l: iLetelier fc 01iveiralll987t Ibemoj 
1989'; Lemos fc Letelier"l993"; Bicak, Lynden-Bell fc Katj'iggg': 'Bicak. Lvnde n-BeU fc Pichonlll993l: Lemos fc Letelieilll994l : 
Lemos fc Letcher 1996 ; Gonz alez fc Espitia 2 003; Garcia fc G o nzalez 2004) . Thin rotating disks that can be considered as a 
source of the Kerr metric were presented bv lBicak fc Ledvinkal il993l). while rotating disk s with heat flow were studied by 
lOonzalez fc Letelierl (|200(i) . Also thin disks with radial tension (^Go nzalez fc Leteliejll99fll^ . magnetic fields ijLetelieil 11993) 
and magnetic and electric fields jKa tz. Bicak fc Lynden-BcU 199ii.) were considered. The nonlinear superp osition of a disk 
and a black hole was flrst obtained bv lLerno s fc Letelierl (1 19 93ll. Perfect fluid disks with halos were studied bv lVogt fc Letelieil 
i2003l) and charged perfect fluid disks bv lvog^^jeteheil |2004) . The stability of som e Gener al Relativistic thin disk models 
using a flrst order perturbation of the energy-mom entum tensor was in vestigated by IUievi~fc Lctclier (2004) . For a survey 
on self gravitating relativistic disks, see for instance I Karas et alJ ll2004) . 

In the works cited above an inverse style method was used to solve the Einstein equations, i. e., the energy-momentum 



* e-mail: danielvt@ifi.unicamp.br 

f e-mail: letelier@ime.unicamp.br; URL: |http://www.ime.unicamp.br/~letelier/| 



2 D. Vogt and P. S. Letelier 



tensor is computed from the metric representing the disk. Another approach to generate disks is by solving the Einstein 
equations given a source (energy-momentum tensor). Essentialy, they are obtained by solving a Riemann-Hilb ert problem 
and are highly nontrivial. This has been used by the Jena group to generate several exac t solutions of disks jKleinl|l997l: 



Neugebauer fc MeinellHooj Ixiein fc R,iciiteilll999l: iKleinlboOll : iFrauendiener fc KleinlEoOll : iKleinlbooi iKleinlboO.Saf 



KleinI 



\2i , 

Recently lOonzalez fc Letelieil l)2004l) constructed models of static relativistic thick disks in various coordinate systems. 
Although the disks have constant thickness, the matter density decreases rapidly with radius and the 2 coordinate, and in 
principle they also may be used to represent both the disk part and the central bulges of galaxies. 

The aim of this work is to consider more realistic three-dimensional models for the gravitational field of galaxies in 
a General Relativistic context. The distribution of matter fills all the space, but the total mass is finite. Essentially we 
formul ate General Relativistic versions in isotropic coor dinates of the potential-density pairs deduced by Miyamoyo and 
Nagai iMivamoto fc Nagailll975l : lNagai fc MivamotolllQT^) and bv lSatolil Jl98d) . Also the influence of galactic rotation on the 
rotation profile of test particles on the galactic plane is estimated using an approximate form of the Kerr metric in isotropic 
coordinates. 

The work is divided as follows. In Sec. |21we present a special form of the isotropic metric and the components of the 
energy-momentum tensor as functions of the metric coefficients, from which the physical properties of the galactic matter 
distributions can be calculated. We also derive expressions for the tangential circular velocity and specific angular momentum 
of test particles in circular motion on the galactic plane. These results are then applied in Sec. |3 and Sec. |3]to construct 
General Relativistic analogues of two Miyamoto-Nagai models and one galactic model proposed by Satoh. As an example of 
application of these models, in Sec.|K|we calculate a few geodesic orbits and compare them with the Newtonian case. In Sec. 
|S|we deduce the first order effects of galactic rotation on the tangential circular velocity of geodesic motion of particles on 
the galactic plane using an approximate form of the Kerr metric in isotropic coordinates, and apply the result on the above 
mentioned models. Finally, we summarize our results in Sec. 



2 EINSTEIN EQUATIONS IN ISOTROPIC COORDINATES 

Let us consider the particular case of axial symmetric space-times whose metric can be writtem in an isotropic form in 
cylindrical coordinates {t,R,z,ip), 

ds^ = e'-f^'^'c^df^ - e^(«'^) (di?2 + dz^ + i?'d(^2) . (1) 

In order to construct a General Relativistic formulation of the known Newtonian models for the mass distribution of galaxies, 
it is convenient to use a particular form of metric Eq. suggested by generalising the form of Schwarzschild's metric in 
isotropic coordinates 

ds'=(^i-^j d'At' + ff {AR" + Az^ + B?d^'') , (2) 
where / = f{R,z). Metric Eq. @ gives the Schwarzschild solution when / is taken as 

(3) 

For metric Eq. ||5J, the Einstein equations G^i, — — ^^T^^ yield the following expressions for the components of the energy- 

= - o , ( + + ^ V (4a) 



t: = ,^aii+mi-f) (^^'- + ^ + 2/,^. - f%) , (4c) 

{.flR.-3fM, (4d) 



47rG(l + /)5(l-/) 



= 4^G(i + /)5(i_j) [/ + /-) - - ] • 

The energy density is given by cr = r*/c^, and the stresses (pressures or tensions) along a particular direction read Pi = — 
when the energy-momentum tensor is diagonal. It is a surprising fact that the component is proporcional to the usual 



^ The Einstein tensor was computed using MAPLE and the result was checked using REDUCE. 
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Laplacian of the function / in flat cylindrical coordinates. Note that in the Newtonian limit when f <^ 1, Eq. Mall reduces 
to Poisson equation 

V^* = 47vGpN, (5) 
if the function / is related to the gravitational potential $ by 

The energy-momentum tensor will be diagonal {T^ — — 0) provided / has the form 

(J 

^ " ^w{R)+g{zy 

where C is a constant and w{R) and g{z) are arbitrary functions. Furthermore, the components and will be equal 
(isotropic radial and azimuthal stresses) only if w{R) = R^. 

Two physical quantities of interest related to the circular motion of test particles along geodesies on the galactic plane 
are the tangential velocity Vc (rotation profile) and the specific angular momentum h. The assumption of geodesic motion is 
valid for the case of a parti cle moving in a very d iluted gas like the gas made of stars that models a galaxy disk. For metric 
Eq. Q the expressions are JVoet fc LeteUeJbooall 

2 2 r,2 e^(e'''),-R /qX 

Vc = c R — , „„ ,, — , (8) 



, _ r2 A / (e'Q.fl , . 

V e-{R2e^).R - i?2e^(e-),fl ' ^ ' 

where all functions are evaluated on 2; = 0. The angular momentum can be used to determine the stability of circular orbits on 
the g alactic p lane by using an extensio n of the Rayleigh criteria of stability of a fluid at rest in a gravitational field jRavleighl 
I1917I see also lLandau fc Lifshitz|[l98i1) 



, dh 
""dR 

For the specific metric form Eq. Eq. simplify to 

2 'ic^Ri.i 



(l-/)(l + / + 2i?/,fl) 



> 0. (10) 

z=0 



(11) 



3 GENERAL RELATIVISTIC MIYAMOTO-NAGAI MODELS 
3.1 First Model 

The simplest gravitational potential proposed by Miyamoto and Nagai jMivamoto fc NagaJllO?^ : iNaeai fc MivamotolllOT^ 

to represent stratifications of mass in the central bulges and in the disk parts of galaxies is given by 

HR, A = — — . (13) 

yi?2 + (a + V^2 + 62)2 

where a,h are positive constants. The corresponding three-dimensional density distribution pn(R,z) is easily derived from 
Poisson equation ((KJ 

Pn[R-,z) = ——- ^ -■ (14) 

[i?2 + (a + + 62)2] (^2 _^ 52)3/2 

Now we consider the function f{R, z) defined in the previous section as 

f{R, z) = =■ (15) 



2c^ J R^ + [a + y/lP^T¥Y 
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Figure 1. Surface plot and level curves of the energy density cr Eq. I16ai and "effective Newtonian density" p Eq. 1171 with parameters 
a = 1 and b = 2. 



Substituting Eq. 11511 in Eq. 14al - ll4ell we find the following expressions for the nonzero components of the energy-momentum 
tensor 



Pr=Pv = 



l + JK^ + {a + iy 



12 [ai?2 + (a + ^)2(a + 2C)' 






l + ^i?2 + (a + C)2 


5 


-1 + ^J + {a + iY 



47rC2 



1 + J + {a + (,y 



l + \/i?2 + (5 + ^)2 



(16a) 



(16b) 



(16c) 



where the variables and parameters were rescaled in terms of the Schwarzschild radius in isotropic spherical coordinates 
Ta = GM/{2c^Y: R = R/rs, z = z/rs, a = a/vs, b = b/rs, a = ^a, p, = ^^pi and ^ = ^/pTIp. The "effective Newtonian 
density" p — a +PrIc^ +P'filc^ Pz/c^ reads 



P = 



b'^R? + (~a + iY 


hR^ + {h + if{a + iC> 






1+ ^J + {a + iY 


5 


-l + ^JR^ + {h + iY 



(17) 



with P= ^p- 

The strong energy condition requires that p > 0, whereas the weak energy condition imposes the condition a > 0. The 
dominant energy condition requires \pR/a\ < (? , \pz/a\ < (? and \pip/a\ < (? . Eq. p6a^ - l|17|l show that the densities are 
positive everywhere and so are the stresses (pressures) if we impose the condition a -\- b > 1. Moreover, we have isotropic 
pressures along the radial and azimuthal directions. When R, z, a, b^ 1, the first order expansion of Eq. 11711 is 



b^ \aR^ + {a + 0\a + 30] 

ri ] 5/2 

47re^ [i?2 + + e)2j 



+(5-1-^)2 



(18) 



Thus the first order General Relativistic correction lowers the mass density. 

In Fig. [Hand Fig. H we show the surfaces and level curves of the densities and pressures Eq. I16at " 117ll for parameters 
a = 1, & = 2 and the same physical quantities for parameters 5 = 1, 6 = 1 are displayed in Fig. As in the Newtonian 

case, the less the ratio b/a, the fiatter are the mass distributions and the pressure distributions show similar behaviour. For 



2 Note that the Schwarzschild radius rs in isotropic spherical coordinates is related to Schwarzschild radius r'g in canonical spherical 
coordinates by = r'^/A. 
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R 

Figure 2. Surface plot and level curves of the radial pressure pjj Eq. <16bt and vertical pressure pz Eq. I16cl with parameters 5 = 1 
and 6 = 2. 




Figure 3. Surface plot and level curves of the energy density a Eq. I16ai and "effective Newtonian density" p Eq. I17i with parameters 
a = 1 and b = 1. 



a = 1, 6 = 2 we have Pr/(j < 0.2 and Pz/o < 0.15 and for 5 = 1, 6 = 1 the ratios are Pr/S < 0.4 and Pz/a < 0.3. Thus for 
these parameters all energy conditions are satisfied. 

The tangential circular velocity Eq. (IIH and specific angular momentum Eq. 11211 for / given by Eq. 11511 evaluated on 
2 = read 



c 



2 



ft= ^ ' ^ ^= , (20) 

(^2 + 772)3/4. / (^2 _^ ^2)2 + ^2 _ ^2 _ 4^2 Y^r^ 



where we defined ft = crgh and r; = a + b. The first order expansion of Eq. 1191 reads 

Vc RV2 L , 



c 



+ (a + 6)2j I ^/J2 _^ (a _^ 5)2j 



(21) 
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Figure 4. Surface plot and level curves of the radial pressure pjj Eq. <16bt and vertical pressure pz Eq. I16cl with parameters 5 = 1 
and 6 = 1. 



Fig.|K|(a)-(b) shows, respectively, the rotation profile Eq. 11911 and the curves of specific angular momentum Eq. 12011 for 
parameters o = 1 and b = 2, b = 1 and b = 0.5. In Fig. 13 (a) the dashed curves represent the newtonian tangential velocities. 
Relativistic tangential velocities are always higher than the Newtonian ones, as suggested by the first order correction in Eq. 
I|21|l . The curves of specific angular momentum show that circular orbits on the plane z — turn to be unstable for higher 
flattened matter distributions (the curve for b = 0.5). The orbits for the parameters a — 1, b — 2 are all stable, whereas in 
the case 6 = 1 a small region of instability appears between R~ 5 and 7? « 6. The Newtonian limit of Eq. l|20|l gives 

-h= ^ (22) 

[7?2 + (a + &)2]3/4' 

from which is easily checked that ^ > for all a, b and R. Thus the instability is a pure relativistic effect ( the same is 
verified for cir cular orbits around the superposition of black holes and axially symmetric structures in Weyl coordinates, see 
iLeteliej toO^ ). 



3.2 Second Model 



Another pair of potential and density functions presented by Miyamoto and Nagai as a generalization of Toomre's model 2 is 
described by 



<E.(i?,z) = 
Pn{R, z) = 



GM 



Vi?2 + (a + C)2 
36^ M 



1 + 



a(a + C) 



i?2 + (a + C)2. 

[R\a' + + (C + afiC + 4< + a)] 



(23a) 



(23b) 



AtiC [R^ + {a + CY]"'' 

where ^ = + 6^. Using / defined by Eq. @ in Eq. 14all - 14ell . we find that the components = are nonzero. The 
physical variables of the matter distribution are obtained by solving the eigenvalue problem for 

T^S'' = AH". (24) 

We find that r"' can be put in the form 

Tab a b , a b , a b , a b /nr\ 

= 0-6(0)6(0) +^+6(1)6(1) +^-6(2)6(2) +Pve(3)e(3), (25) 

where 

'! + / 



= (0) 



P± 



Tl + Tl 
2 



1-/ 



0,0,0 , 



+ 4(r/ 



°W = (0,e(i),e(i),0), ef2) = (0, 6(2) , 6(2), 0), 

1 



_rplf 



H3) 



°'°'%(1+/)^ 



(26) 



(27) 
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and 



(i + /)^vWFTTTfT^' ''^ " (i + /)V(W+TTf+p^' 

The effective Newtonian density reads p = a + p+ / (? + p- / (? + / (? — a — T^/c^ — /c^ +p^/c'^ ■ The explicit expressions 
for a, p and p^p are 

36^X' iR^a' + e) + (a + O'ia' + + 4501 
= Ts (29) 



47r53 + R2 + (a + C)(2a + C)] ' 

[x^/^ + + (a + e)(2a + o] ' [x^/^ „ 7?^ - (a + 0(2a + 0] ' 



s&^x"/' [^'(a' + C') + (a + C)'(a' + + 4ac)l 

P= -r 7' (30) 



^4 . r 



Stt^s [x3/2 + + (a + ^) (25 + ^)] [^3/2 _ 7^2 _ (a + ^) (25 + ^)] 
+R^{d. + [ga^fc^ + UaP^"^ + 4pf + 25^(76^ + + {a + 0"* [6a^&^ + 1236^^'^ + 26^^^ + a^(19&^ + } , (31) 
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Figure 6. Surface plot and level curves of the energy density a Eq. 1291 and "effective Newtonian density" p Eq. 1301 with parameters 
a = 1 and b = 2. 
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Figure 7. Surface plot and level curves of the pressures p± Eq. I26i with parameters a = 1 and 6 = 2. 



where ^ — \/z^+lP, x = + {o. + 0^ '•^^'^ the dimensionless variables were previously defined. The components T^, 
and read 



SttC^ + f2 + (a + ^)(2a + e)]' [x'/2 „ ^2 _ (a + ^)(2a + ^)] 



+f^{a + [ga^P i25.Pf + iP^^ + 2h^(i? - 52^)^] + (a + if [ea^P + i2ai?i^ + 2p^^ + a^(i9p + 5^)^] | , (32) 



rl 47re [x'/2 + f2 (a e)(2a + C)]' [x^^' - f 2 _ + ^)(2a 0] ' 



(33) 



rl + ^2 (a _^ ^)(2a + [^3/2 _ ^=2 _ (g + ^)(2a + c)] 

{^^^'(a' - - 2r'?(a + 0' [356^ + C(2a' + b')] + (a + 0^ - 86") - b^i{6a + } , (34) 

and p± is further adimensionalized as p± — ^^^p± ■ 

The condition (a + 6)^ > (2a + b) ensures positive densities and non-singular behaviour for the densities and pressures. 
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Figure 9. (a) Relativistic rotation profile Eq. 1361 and Newtonian rotation profile (dashed curves) with parameters a = 1 and b = 3, 
6 = 2. (b) Specific angular momentum Eq. 1371 for the same parameters. 



The first order expansion of Eq. lOnj reads 

_^ 3P[R^(a^ + e) + {d + OHa'+e+4~aO] ( 4[fi2 ^ (g + e)(2a + Q] ] 
47rC3[i?2 + (a + 02]''/2 I [i?2 + (a + 0']^/' i' 

Tlie surfaces and level curves of tlie densities Eq. I|29|l - H3fl|l and pressures Eq. I)26|l are siiown in Fig.|^and Fig.Q respectively, 
for parameters a = 1, b — 2. Fig. |HJa) displays the surface and level curves for the azimuthal pressure. For these parameters 
we havep,p/cr < 0.3, p+/a < 0.25 and p-ja < 0.3. 
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The tangential circular velocity Eq. IIH and specific angular momentum Eq. 1121 are given by 

(36) 







[x»/2 - ^2 - 77(25 + b) 





V2R^ \x'/' + R'+ v{2a + b)] Jr^+ ri{4a + b) . . . ^ ,^ . 

= ^ ^9/4 ^ |x + [i? + '7(25 + b)\ [r" + 5ar,R^ - 7?^(2a + b)\ 

-AR\'/^[R^+vi^a + b)\y'^\ (37) 



where h — crsh, r] — a + b and x ~ R^ + ■ The first order expansion of Eq. 1361 reads 

,1/2 



^27? [i?2 + (a + 6)(4a + &)] j ^j?2_^(a + 6)(4a + 6) 



5/4 1 1 + —JZ _ .5/2 ^ ^ • (38) 



[j?2 + (a + 6)2j 1^ Ji^2 + (a + ^)2j 



In Fig. |nia)-(b) we display the rotation profiles and curves of specific angular momentum for parameters a = 1, b = 3 and 
b = 2. As in the previous case, relativistic effects increase the tangential velocity. We also note that the orbits for b = 3 and 
b = 2 are all stable. 



4 A GENERAL RELATIVISTIC SATOH MODEL 

ISatohl l)l98(tl derived other pairs of three-dimensional potential and density functions by inflating the Plummer models of 
order n. In the limit when n ^ oo the pair takes the simple form 

Gm 



<i>{R,z) = 
Pn{R,z) = 



i?2 + ^2 _^ „(a + 2Vz^ + b2) 

b^AI [R^ + + (g + 2Vz2 + b2)(a + 3V22 4. 52)] 



(39a) 
(39b) 



4^ (22 + 62)3/2 [^2 j^^2j^ a{a + 2V22 + b2)] '''' ' 

Using the function f(R, z) with <E> given by Eq. ll^ . we get the following expressions for the nonzero components of the 
energy-momentum tensor Eq. I4al - ll4ell 

aP + + + 2^(5 + 30 



PR =P^ 



l + jR^ + P+aia + 2^) 









l-f ^^^2 + 52 + ^(5 + 2^) 


5 


-l + ^i?2 + 52 + a(a + 20 



aP[a + 2i) 



47r^2 



1 + J R^ + P +h{a + 2i) 



1+ J R^ + + a{a + 2£,) 



P = 



aP ^ + + h{a + 2i) 


+ 52 + (5 + 2^(5 + 30 




47rC^ 


1 + ^ + +h{a + 2£) 


1 


-i + \/7?2 + 52 + a(a + 20 



(40a) 



(40b) 



(40c) 



(40d) 



where ^ — \J + }p- and the dimensionless variables were previously defined. As / has the form of Eq. Q, the energy- 
momentum tensor is diagonal and the radial and azimuthal stresses are equal. We have ct > and p > everywhere and 
pressures without singularities provided 5(2 + 2b) > 1. In the weak field limit, Eq. (I40dt may be expanded as 



ab 



2 [i?2 + 52 + (5 + 20(5 + 30] 



47r 



e \r^ + 52 + 5(5 + 20] 



5/2 



fl;2 + 52 + a(a + 20 



(41) 



which introduces again a negative first order relativistic correction to the Newtonian mass density. 
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Expressions for the tangential velocity and angular momentum follow directly from Eq. lllll - 112t 



h= , ^ ^ (43) 



(44) 



where /i = crs/i and r] = S(5 + 26). The first order expansion of Eq. 14211 is given by 

RV2 , R^ 



C r - - 1 3/4 1 ' r . - 1 3/2 

i?2 + fi(a + 26) y i?2 + a(a + 26) 



The surfaces and level curves of the densities and pressures Eq. 14Uall - ll40dl l for parameters 2 = 1, 6 = 2 are plotted in 
Fig. 1101 and Fig. 1111 We have pr/(t < 0.6 and Pz/o" < 0.5, so all energy conditions are satisfied. Rotation profiles and curves 
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of specific angular momentum are displayed in Fig. If 21 for parameters d — 1, b — 2 and b — 1. We note that circular orbits 
are stable for b = 2. 



5 GEODESIC ORBITS 

An interesting application of the General Relativistic galactic models is the study of geodesic orbits outside the galactic plane. 
As an example we calculate numerically two orbits for the first Miyamoto-Nagai model (Sec. 13. ft and compare them with the 
Newtonian case. 

We solve numerically the geodesic equations of motion 

+ r^^i^i" = 0, (45) 

for metric Eq. ||5J, where FJ^^ are the Christoffel symbols and the dot denote differentiation with respect to the proper time. 
With the orthonormal tetrad defined as 

^+1(1,0,0,0), eln-) = j^^j^{0,l,0,0), 



1 - 

f „ , , f 



(46) 



^(^) = (TTTF^"' °' ^' " j^(i + /)2 (°' 0' 0' 1)' 

the tetrad components of the four-velocity read 

= 7(f , V sin ij) cos 6, v sin sin 9, v cos tp), (47) 



with 7 = — (we have set c = f). The specific energy and angular momentum of the test particle are 

f = 7^, (48) 
h = r{l + ff-yvcosij. (49) 
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(a) 15 




Figure 13. (a) Geodesic orbits for the first Miyamoto-Nagai model with parameters a = 1, fe = 1. Initial conditions: Rq = 15, ^p = 45° 
with energy £ 0.901 and angular momentum h ~ 5.03. (b) Newtonian orbit with energy £ — I and same angular momentum. Initial 
conditions: Ro ~ 19-84 and V ~ 36.37°. 




012345 012345 
ip/2jl 9/2JI 



Figure 14. The coordinates R and z as functions of the number of revolutions ip/2TT for the orbits shown in Fig. 1131 (a) Geodesic and 
(b) Newtonian orbit. 

As initial conditions we take a position at radius _Ro on the disk's plane and components of the four- velocity = 7(1, 0, vq sin ip, vg cos ip) , 
where vq is equal to the tangential velocity of circular orbits at radius Rq. Similar initial conditions for the Newtonian case 
are chosen with the same angular momentum and energy £ — 1. Both the geodesic and Newtonian equations are solved with 
the fourth-fifth Runge-Kutta method with adaptive size step. 

In Fig. ll3ttT^ we compare both orbits for the first Miyamoto-Nagai model with parameters a = 1 and b — 1. For the 
orbits shown in Fig. llStlTll large radii have been chosen as initial conditions. As expected, General Relativistic effects are 
small. A noticeable difference is seen in a periodic variation of the amplitude of the R coordinate as function of the azimuthal 
angle in the relativistic case. On the other hand, the orbits shown in Fig. 11514161 have been computed near the center where 
the gravitational fields are strong, and the deviation between relativistic and Newtonian orbits are much more pronounced. 

6 FIRST ORDER EFFECTS OF GALACTIC ROTATION ON THE ROTATION PROFILES 

The galactic models presented in the last sections are all static. Astrophysical objects are expected to have angular momentum, 
thus more realistic models should incorporate rotation. An exact stationary General Relativistic galactic model could be 
constructed, for example, if one managed to find a closed form of the Kerr metric in terms of elementary functions in 
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(a) 0.8 




Figure 15. (a) Geodesic orbits for the first Miyamoto-Nagai model with parameters a = 1, 6 = 1. Initial conditions: Rq = 1, ip = 45° 
with energy £ ^ 0.196 and angular momentum h Ri 0.867. (b) Newtonian orbit with energy £ — 1 and same angular momentum. Initial 
conditions: Ro » 1.88 and tp 38.01°. 
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Figure 16. The coordinates R and z as functions of the number of revolutions ip/2iT for the orbits shown in Fig. 1151 (a) Geodesic and 
(b) Newtonian orbit. 



stationary isotropic coordinates, and then repeat the procedure outlined in the previous sections to calculate the matter 
properties. Although this cannot be done for the exact Kerr solution, the same is not true for the first order expansion of the 
Kerr metric with respect to the rotation parameter. It reads 

2GM\ 2, ,2 dr^ 2.,.2 , . 2., 2^ , 4GMa 



{de^ + sin^ edip^) + — — - sin Odtdip, (50) 



rc^ / 1 — 2£M 

where a is the rotation parameter. If we transform metric Eq. 1501 to cylindrical isotropic coordinates 
(i, R, z, ip) and express it in terms of the function / given by Eq. igj we get 

e^di^ - (1 + {dR' + dz^ + R'd^'') + ^^^^dtd^. (51) 

Now the function / may be replaced by any of the forms presented in the previous sections to generate spacetimes with 
matter. However, a direct calculation using Einstein's field equations and metric Eq. 1)5 1^ shows that the influence of rotation 
on the physical variables of the matter is of at least 0{a^), and thus to have a consistent correction one should also start 
with metric form Eq. H51|l corrected at least to 0{a^), but it does not seem possible to express this particular form of the 
metric with second order corrections of the rotation parameter. Despite this, a consistent first order correction of the rotation 
profiles can be calculated. 
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We rewrite Eq. 1511 as 

dr^ = gttdt'^ - Qrr [AR" + Az^ + R^Aip^) + 2gt^AtAip, (52) 

where for convenience we liave set c = 1. Let t;" = {t,0,Q,tp) be the velocity four- vector, where dots denote differentiation 
with respect to r. The geodesic equation for circular orbits is 

t ^ 



gtt.Rt + '2gtip,Riip ~ {R qrr), Rip 



0^¥ = 



[R^gRR),R 

Thus the prograde and retrograde angular velocities are in general different. We define an orthonormal tetrad e[p)i as 



(53) 



/ 5?^ + R^gttgRR 



,0,0,0 



Ry/gRR 

e(R)t = (0, -^JgRR, 0, 0), e(2)i = (0, 0, 



/ffifl, 0), 



gtf 



, R^/gRR 

and project the components of on the tetrad: v'"' = 



,0,0, -Ry/gRR 



(a) i 

e) 'v 



(t) 



gL + R^gttgRR 



RJWrr 



gtv 



Using Eq. (1^ in Eq. (I55al " ll55bll we get 

- 1 

v(0 ^ 



3*1= + R'^gttgRR 



Ry/gRR 

R^gRR 



i + R^/gRR(p. 



gtf + ,r}0 {gtv,R ± \/ {gtv,RY + gtt,R{R'^gRR).R) 

(H^gRR).R \ I 



(54) 

(55a) 
(55b) 

(56) 



The terms g1^ and ijjt^.R)^ in Eq. I56II are of 0{o?^ and thus can be neglected. Using metric Eq. 1511 . Eq. 1561 can be 
rearranged as 



(1 -/)(! + , f + 2fl/,fl) 



±1 + . 



2^afRc' 



-Rf.Ri^ - f) 



(57) 



G2A/2(l + /)3 y 2(1 + / + 2i?/,H) 

where c has been reintroduced. The term outside brackets is the tangential velocity Eq. 1111 . The term with a is always 
positive, thus the first order effect of rotation increases (decreases) the prograde (retrograde) tangential velocity. 
Eq. 1)57^ calculated with the functions / used in Sec. Em results m 



Vc 

c 
c 



V2RiR^ + 0^ 



/4 



-1 + (i?2 + ^)1/2]1/2[_7J2 + ^ + (iJ2 + ^)3/2]l/2 



±1 + 



^/2i?C^/''[i^2 + {a + b){4a + b)^ 



3Qfi2(i?^ + C) 

[i + (i?2 + ^)i/2 

/2 



1/4 



2[-l + (fl2+g)l/2] 
-i?2 + ^ + (i?2 + ^)3/2 



(58) 



1^(^3/2 _ 7^2 _ (5 _^ 6)(2a + 6)] [C'/^ - - 55(5 + 6)7?2 + (a + &)3(2o + 6) 



1/2 • 



±1 + ' 



■i&if(*/*[R^ + (a + 6)(2a + b)f 2[R2 + (a + b)(4a + 6)][C3/2 _ 7?2 _ _^ ^)(2a + &)] 



(59) 



[^3/2 + R2-^-(^a + b){2a + b)]'^ V c^''^ - i?4 _ 5a(a + 6)i?2 + (a + fo)3 (2a + b) 

where in Eq. 1581 ^ = (a + 6)2 for the first Miyamoto-Nagai model (Sec. l3.ll . (, — a{a + 2b) for the Satoh model (SecllJ, and 
in Eq. 1591 = n + {a + b)^ for the second Miyamoto-Nagai model (Sec. 13.21 . The dimensionless variables were defined in 
Sec. 13. H and a = a/vs- 

In Fig. I17f a')-(b') we plot the curves of progade and retrogade tangential velocity for the first Miyamoto-Nagai model 
with parameters a = 1, & = 2 and a = 1, fe = 1 with rotation parameter a — 0, a = 0.2 and a = 0.4. The same curves for the 
second Miyamoto-Nagai model and for the Satoh model with parameters a — 1, b = 2 are shown in Fig. ^Ja)-(b). 



7 DISCUSSION 

Using a special form for the metric in cylindrical isotropic coordinates, we have constructed what may be interpreted as the 
General Relativistic versions of some known Newtonian models for th r ee-dimensional distributio ns of matter in gal axies, like 
those formulated by Miyamoto and Nagai jMivamoto fc Nagailll97,4 iNagai fc Mivamotolll976^ and ISatohl l)l98d) . For the 
first Miyamoto- Nagai potential and for one Satoh potential we obtained matter distributions that have similar properties: 
their energy-momentum tensor is diagonal with isotropy of pressure in the radial and azimuthal directions; energy density 
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and the "effective Newtonian density" are positive everywhere as well as the stresses (pressures) in the vertical direction. For 
particular values of the free parameters we found that the dominant energy condition is also satisfied. Stability analysis of 
the circular geodesic orbits on the galactic plane using an extension of the Rayleigh criteria of stability shows stable orbits 
for some combination of free parameters. The second Miyamoto-Nagai potential function, after diagonalization of the energy- 
momentum tensor, yield matter distributions with positive energy density and three different pressures. Also in this case, free 
parameters can be chosen so that all energy conditions are satisfied and circular orbits are stable on the galactic plane. 

As an example of application of the models, we have numerically calculated some geodesic orbits for one of the potentials 
and compared them with the Newtonian orbits with the same energy and angular momentum. Near the central regions where 
the gravitational fields are strong, the motion of particles is considerably altered by General Relativistic effects. 

We also calculated the first order effects of galactic rotation on the tangential velocity of circular orbits on the galactic 
plane using an approximate form of the Kerr metric expressed in cylindrical isotropic coordinates. In general, rotation increases 
the progade tangential velocity and has an opposite effect on the retrogade tangential velocity. 

It should be mentioned that the stability study of the models presented based on the extension of Rayleigh criteria of 
stability is very limited. A more realistic stability analysis should rely on the first order perturbed Genera l Relativistic fluid 
equati ons taking into account two spatial coordinates, which may be not a trivial task (see, for instance, lUievic fc Letelieil 
J200J) for the one dimensional case). This will be subject of further investigation. 
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Figure 18. Prograde (upper) curves and retrogade (lower) curves of tangential velocity for (a) the second Miyamoto-Nagai model Eq. 
I59i with parameters a = 1, 6 = 2 and (b) for the Satoh model Eq. I58i with parameters a = 1, b = 2. Rotation parameter: a = (solid 
line), a = 0.2 and a = 0.4. 
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